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Resilience, a system’s ability to adjust its activity to retain its basic 
functionality when errors, failures and environmental changes 
occur, is a defining property of many complex systems 1 . Despite 
widespread consequences for human health 2 , the economy 3 and 
the environment 4 , events leading to loss of resilience — from 
cascading failures in technological systems 5 to mass extinctions 
in ecological networks 6 — are rarely predictable and are often 
irreversible. These limitations are rooted in a theoretical gap: the 
current analytical framework of resilience is designed to treat 
low-dimensional models with a few interacting components 7 , 
and is unsuitable for multi-dimensional systems consisting of 
a large number of components that interact through a complex 
network. Here we bridge this theoretical gap by developing a set of 
analytical tools with which to identify the natural control and state 
parameters of a multi-dimensional complex system, helping us 
derive effective one-dimensional dynamics that accurately predict 
the system’s resilience. The proposed analytical framework allows 
us systematically to separate the roles of the system’s dynamics 
and topology, collapsing the behaviour of different networks 
onto a single universal resilience function. The analytical results 
unveil the network characteristics that can enhance or diminish 
resilience, offering ways to prevent the collapse of ecological, 
biological or economic systems, and guiding the design of 
technological systems resilient to both internal failures and 
environmental changes. 

The traditional mathematical treatment of resilience used 
from ecology 4 to engineering 5 approximates the behaviour of a 
complex system with a one-dimensional (ID) nonlinear dynamic 
equation 

£=/gm (i) 

at 

The functional form of f(/3, x) represents the systems dynamics, and 
the parameter (3 captures the changing environmental conditions. The 
system is assumed to be in one of the stable fixed points, x 0 , of equation 
(1), extracted from 5 

f(/3,x 0 ) = 0 (2) 


where equation (2) provides the systems steady state and equation (3) 
guarantees its linear stability. The solution of equations (2) and (3) pro- 
vides the resilience function x(/3 ), which represents the possible states 
of the system as a function of / 3 (Fig. la-c). The shape of this function 
is uniquely determined by the functional form of f(/3, x). In contrast, 
the momentary state of the system is determined by the tunable param- 
eter f3. At some critical point (3 C the resilience function may feature a 
bifurcation (Fig. la) or become non-analytic (Fig. lb, c), indicating that 


the system loses its resilience by undergoing a sudden transition to a 
different 8,9 , often undesirable, fixed point of equation (1). 

Although it is conceptually powerful, this analytic framework does 
not account for the exceptionally large number of variables that in 
reality control the state of a complex system. Indeed, real systems 
are composed of numerous components linked via a complex set of 
weighted, often directed, interactions 10,11 , and controlled by not one 
microscopic parameter, but by a large family of parameters, such 
as the weights of all interactions. Hence, instead of a ID function 
f((3, x), characterized by a single parameter f3 y their state should be 
described by a network of coupled nonlinear equations that cap- 
ture the interactions between the system’s many components, and 
account for the complex interplay between the system’s dynamics 
and changes in the underlying network topology 6,12 . The resulting 
resilience function is therefore a multi-dimensional manifold over 
the complex parameter space characterizing the system (Fig. ld-f), 
which, using the current tools, cannot be treated analytically. Here 
we overcome these longstanding limitations by developing a general 
network-based theoretical framework that allows us to explore and 
predict the multiple roots and dimensions of resilience, exposing cru- 
cial determinants of resilience loss in complex natural and man-made 
systems. 

Consider a system consisting of N components (nodes) whose activ- 
ities x = (x b . . . , x n Y follow the coupled nonlinear equations 12,13 

dx- N 

-± = F(x i ) + J2A ii G(x i ,x j ) (4) 

d t j - 1 

The first term on the right-hand side of equation (4) describes the 
self- dynamics of each component, while the second term describes 
the interactions between component i and its interacting partners. The 
nonlinear functions E(x z ) and G(x z , xj) represent the dynamical laws 
that govern the system’s components, while the weighted connectiv- 
ity matrix Ay captures the interactions between the nodes. With an 
appropriate choice of F(x z ) and G(x z , x ; ), equation (4) is used to model 
numerous systems known for their resilience, ranging from cellular 14 
to ecological 15,16 and social systems 17 . 

In analogy with the ID system of equation (1), a transition from a 
desired to an undesired stable fixed point captures the loss of resilience 
in the multi- dimensional system of equation (4). The key difference is 
that in equation (4) resilience loss can be induced by changes in any of 
the N 2 parameters of the weighted network Ay, each change capturing 
a different kind of perturbation (Fig. lg). For instance, the extinction/ 
introduction of species in an ecological system may correspond to the 
removal/addition of one or several nodes 7,18 . Alternatively, the loss of 
an enzyme catalysing a reaction in a metabolic network 19 may corre- 
spond to link removal. Finally, global environmental changes, such 
as increases in ocean acidity or global warming 20 , may correspond to 
global shifts in the weights of Ay. 

We illustrate the challenges such multi-dimensional systems offer by 
exploring the mutualistic interactions among species in an ecological 
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Figure 1 | Network resilience, a-c, In ID systems 
resilience is captured by the resilience function 
x([3), which describes the state(s) of the system as a 
function of the tunable parameter (3. We illustrate 
three typical examples, a, The bifurcating resilience 
function. The system exhibits a single stable fixed 
point for (3> (3 C (blue) and two (or more) stable 
fixed points, a desired (blue) and an undesired 
(red) for /3 < (3 C . b, Resilience function with a first- 
order transition from the desired (blue) state to the 
undesired (red) state, c, Resilience function with a 
stable solution for (3 < (3 c and no solution above (3 C , 
resulting in an uncontrolled divergent or chaotic 
behaviour, d-f, In a multi- dimensional system, the 
single parameter f3 is replaced by the complex 
weighted network Ay, whose characteristics 
depend on both environmental conditions and 
the specific pairwise interaction strengths. 
Consequently, the resilience function, now 
capturing the behaviour of the vector state x(Ay), 
is a multi- dimensional manifold, prohibiting 
analytical treatment. The 3D plots show the 
resilience plane for a four- node system, displaying 
x(A 2 3, A24) for fixed A 12 and A 34 . The full 
description of an IV-dimensional system requires 
an N 2 - dimensional plane, tracking the state of the 
system as a function of all Ay. g, After applying our 
formalism, the multi- dimensional manifold shown 
in d-f collapses into a ID resilience function in 
(3 - space (blue and red solid lines). The structure of 
this function, and hence its critical points (3 c ef{ 
(dashed lines) is fully determined by the systems 
dynamics F(xj) and G(x iy xj): ecological, regulatory, 
power transmission and so on (left). The network 
topology Ay (right) determines /3 e ff through 
equation (8), and hence the specific state of the 
system (brown dot) along the resilience function. 


network. Here equation (4) tracks the abundance Xj(t) of species /, 
following 16 


i£=j, + xi 

d t 


1 - — 

K, 


* -1 




N 

+ 2 ^ij 

j = i Dj + Ej Xj + Hj Xj 


( 5 ) 


The first term on the right hand side of equation (5) accounts for the 
incoming migration of i at a rate Bj from neighbouring ecosystems. 
The second term describes logistic growth with the system carrying 
capacity 21 K b and the Allee effect, according to which for low abun- 
dance (Xj < Q) the system features negative growth 22 . The third term 
describes mutualistic interactions, captured by a response function 
that saturates for large Xj or Xp indicating that / s positive contribution 
to Xj is bounded. To construct Ay we used symbiotic interactions, 
such as plant-pollinator relationships, collected for seven ecologi- 
cal systems 23 (Supplementary Table 1), describing networks ranging 
from N= 10 to N = 1,429 nodes (Fig. 2a, b). We numerically solved 
equation (5), and tested its resilience under three realistic pertur- 
bations (Fig. 2c): first, we randomly removed a fraction f n of nodes, 
capturing plant extinctions; second, we removed a random frac- 
tion f\ of pollinators, thus perturbing the mutualistic link weights; 
and finally, we randomly rescaled all the weights Ay, reducing their 
strength on average by a fraction/ w to mimic global environmental 
changes. 

We find that for small perturbations the system maintains its resil- 
ience: its only stable fixed point is * H , in which the average abun- 
dance (x) is high. However, when the perturbation exceeds a certain 
threshold a bifurcation occurs, resulting in two stable fixed points: the 
desired x H and an undesired low-abundance state x L (Fig. 2d). Under 
these conditions the system loses its resilience, potentially transi- 
tioning to the undesired x L . The precise bifurcation point marking 
this loss of resilience is, however, highly unpredictable. For instance, 


Netl displays a different resilience pattern for node removal (Fig. 2d), 
link removal (Fig. 2e) or global weight changes (Fig. 2f), indicating that 
different forms of perturbations lead to different outcomes within the 
same system. Such diversity is also observed for similar perturbations 
across different systems. For example, while we need to remove at least 
35% of the pollinators for Netl to lose its resilience (Fig. 2e), Net5 
remains resilient even after 80% of its pollinators are deleted (Fig. 2h). 
Finally, even the sequence in which the perturbation is applied makes 
an important difference: Netl can lose its resilience anywhere between 
the removal of 30% to 70% of its nodes (Fig. 2d), depending on the 
specific trial. 

Altogether, we analysed fourteen experimentally mapped networks 
(Supplementary Table 1), finding that the resilience transition depends 
on the network topology, the form and the nature of the perturbation 
applied and the specific realization (Fig. 2d-l, and Supplementary 
Figs 2-4), exposing severe limits to our ability to predict the network 
resilience. We now show that this seemingly unpredictable behaviour 
can be systematically treated by focusing on the systems natural state 
and control variables. In a network environment, the state of each 
node is affected by the state of its immediate neighbours. Therefore, 
we characterize the effective state of the system using the average 
nearest-neighbour activity 


^•eff 


1 T Ax 

1 T A1 



(6) 


where s out = (s° ut , . . . , s™ 1 ) 7 is the vector of outgoing weighted degrees 
and s in = (sl n , . - - , 5 iv) T is the vector of incoming weighted degrees, 
(s out x) = (l/N)E|IiS° ut *z, ( 5 ) = (s in ) = (s out ) is the average weighted 
degree, and 1 is the unit vector 1 = (1,...,1) T . As we show in 
Supplementary Information section I, if Ay has not much degree 
correlation, the variable x e ff in equation (6) allows us to reduce the 
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Figure 2 | Resilience in ecological networks. 

a, The bipartite network M& consists of nodes 
representing symbiotically connected species, 
such as plants and pollinators or fish and anemone. 

b, From M& we construct two mutualistic networks 
(Supplementary Information section II.B) by linking 
pairs of plants that share mutual pollinators (Ay), or 
pollinators that share mutual plants (By), c, We tested 
the resilience of fourteen mutualistic networks 
(equation (5); B,- = B = 0.1, Q=C= l,K i = K=5 J 

D t = D = 5, Ei = E = 0.9 and Hj = H = 0. 1 ) against: 

(1) extinction of a fraction f n of plants; (2) extinction 
of a fraction fi of pollinators, impacting the relevant 
plant link weights; (3) decreasing all weights on 
average to a fraction f w of their original value, 
simulating a global change in the environmental 
conditions, for example, varying temperature, d, The 
average abundance of fish in Netl versus f n across 100 
realizations (we highlight one of these realizations in 
black). At a critical fraction f n the system undergoes a 
bifurcation, where in addition to the high- abundance 
state (x H ) an undesired low-abundance state emerges 
(x L ). However, owing to the multi-dimensionality of 
Ay, each realization is microscopically distinct, and as 
a result the bifurcation point is unpredictable, ranging 
from/ n = 0.3 to f n = 0.7 across different realizations, 
e, f, Similar diversity characterizes the systems 
response to link perturbation/i (e) and global 
perturbations/ w . (f). g— 1, The difficulty in 
predictability is also observed in the larger networks 
Net5 and Net7. Additional results appear in 
Supplementary Figs 2-4. m, Our formalism predicts 
that in /5-space the resilience function takes a 
universal form, similar to that shown in Fig. la, with 
bifurcation at /3 c eS = 6.97, independent of Ay. n, All 
data in d-1 and in Supplementary Figs 2-4, 
comprising 28 highly diverse networks, collapses 
onto the universal resilience function predicted in m, 
indicating that regardless of the network structure and 
the form of perturbation, the state of the system is 
fully determined by (3 e $ (equation (8)). 


multi-dimensional equation (4) to an effective ID equation, written in 
terms of x e ff 

^ = F(x eff ) + /? ef f G(x eff , X eff ) (7) 

dt 

where 


Peft = 


l T As in 


1 t A1 


(s out s in ) 

(s) 


( 8 ) 


averages over the product of the outgoing and incoming degrees 
of all nodes. This reduction maps the multi- dimensional complex 
system (4) into an effective ID equation of the form of equation (1), 
where 


f (/^eff > -^eff ) — F (-^eff ) T fteffG (x e ff , X e ff ) (9) 


of the perturbation (node/link removal, weight reduction or any 
combination thereof), its impact on the state of the system is fully 
accounted for by the corresponding changes in /5 e ff- This implies that 
the rather diverse and unpredictable behaviours observed in Fig. 2 
are, in fact, drawn from a single universal resilience function, inde- 
pendent of the network topology Ay, and uniquely determined by the 
systems dynamics F(x z ) and G(x z , xj). The network Ay, which is fully 
condensed into the single macroscopic parameter /5 e ff> determines 
only the specific state of the system along this function. Such mapping 
of equation (4) to the ID equation (7) allows us to take advantage 
of the theoretical tools developed for low- dimensional systems and 
apply them to a broad range of complex systems. 

To illustrate the power of our formalism we apply it to the mutualistic 
networks of Fig. 2. Reducing the multi- dimensional equation (5) to 
the form of equation (9) we arrive at the ID equation (Supplementary 
Information section II) 


We predict, therefore, that the N 2 parameters of the microscopic 
description Ay collapse into a single macroscopic resilience parameter 
f3 e ff (see equation (8)), so that regardless of the microscopic details 


dv e ff 

dt 


— B + x e ff 


2 ^eff 

~K~ 


(M _ i) + /3 ff *k 

(C ) D+(£ + H)x eff 


( 10 ) 
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Figure 3 | Resilience in gene regulatory networks. 

We ran Michaelis-Menten dynamics (equation (11)) on 
the transcription regulatory networks of S. cerevisiae 24 
and E. coli 25 (Supplementary Table 2) to model the 
dynamics of genetic regulation, providing the average 
activity (x) of all genes, a, By removing a sufficiently 
high fraction f n of nodes the system undergoes a 
transition from x H > 0 (alive) to x L = 0 (cell death). 

The transition point occurs at a different value of f n in 
each realization, b, c, Similar results are found for link 
perturbations and global weight changes, d-f, The same 
diversity is observed for the E. coli network as well, g, Our 
formalism predicts that the behaviour of gene regulation 
is captured by a universal resilience function of the form 
of Fig. lb, with a single first-order transition at /3^ ff = 2. 
This function is determined by the regulatory dynamics 
of equation (11), and is independent of the network 
structure or the nature of its perturbations, h, The results 
of a-f shown in /3-space. Regardless of the systems 
microscopic details, all observed data points, taken 
from S. cerevisiae or E. coli, induced by node 
perturbations (red), link perturbations (green) or weight 
changes (blue), collapse onto the same curve, well 
approximated by the analytically predicted universal 
resilience function (solid line). 


Using equations (2) and (3), equation (10) predicts a bifurcating resil- 
ience function of the form of Fig. la, and a transition from a resilient 
state with a single stable fixed point, x H , to a non-resilient state in which 
both x H and x L are stable. The critical point of this bifurcation is pre- 
dicted to be /3g ff = 6.97, a value fully determined by the dynamics, inde- 
pendent of the network topology Ay (Fig. 2m and Supplementary 
Information section II). 

Our formalism predicts that the diversity observed in Fig. 2d— 1 and 
in Supplementary Figs 2-4, is, in fact, driven by the universal curve 
of Fig. 2m. This universality can be exposed by transitioning to the 
natural parameter space of x e ff (see equation (6)) and /3 e & (see equation 
(8)). Hence we re-plotted all the data of Fig. 2d— 1 (and Supplementary 
Figs 2-4, 28 diverse networks in total), in this effective /3-space, finding 
that, as predicted, all data points collapse into a single universal curve, 
regardless of the size and the topology of the specific ecological net- 
work or the nature of the applied perturbation (Fig. 2n). This collapse 
indicates that our analytically predicted resilience function exposes a 
universality sustained across networks of different size, density, degree 
and weight distributions. Additional extensive testing of this universal- 
ity appears in Supplementary Information section V. 

In summary the resilience pattern of a complex system is effectively 
unpredictable in the natural (x. Ay) state parameter space. Once, how- 
ever, we map the system into /3- space we can accurately predict the 
systems response to diverse perturbations and correctly identify the 
critical points where the system loses its resilience. 

Next we explore the resilience of gene regulatory networks governed 
by the Michaelis-Menten equation 14 


the cell undergoes a transition from a resilient state ((x) > 0) to cell 
death ((x) = 0, Fig. 3a-f). Once again, resilience loss strongly depends 
on the nature of the perturbation (gene knockout, suppression of 
regulatory interactions, environmental change), as well as on the dif- 
ferences between the wiring diagrams of S. cerevisiae versus E. coli 
(Supplementary Table 2). Rewriting equation (1 1) in the reduced form 
of equation (7), we find that regulatory dynamics follow the universal 
resilience function shown in Fig. lb, featuring a single first-order tran- 
sition from the active state to cell death at (Supplementary Information 
section III) 

J 


Indeed, in /3-space, all trajectories of Fig. 3a-f collapse onto the ana- 
lytically derived resilience function, as predicted by our formalism 
(Fig. 3h). In Supplementary Information section IV we develop another 
application area, demonstrating how to apply our formalism to predict 
the resilience of energy supply in the power grid. 

Although the resilience function is uniquely determined by the 
dynamical functions F(xi) and G(x if xj), the actual position of the 
system along this curve, capturing its momentary state, is deter- 
mined by the weighted network topology Ay, as expressed through 
/3 e ff in equation (8). This prompts us to ask what aspects of the 
network topology determine a systems resilience. We therefore rewrite 
/3 e ffas 


d Xi 

d t 


-Bx{+J2A 


'V+i 


(ii) 


The first term on the right-hand side of equation (11) describes degra- 
dation (f= 1) or dimerization (f= 2); the second term captures genetic 
activation, where the Hill coefficient h describes the level of cooper- 
ation in gene regulation 14 . Applying this model with B = 1 ,/= 1 and 
h— 1 to the transcription networks of Saccharomyces cerevisiae 24 and 
Escherichia coli 25 , we find that under sufficiently large perturbations 


Aff — i s ) + Shi ( 13 ) 

where (s), S and 7 ~L represent three characteristics of Ay. The depend- 
ence on the network density (s) indicates that denser networks have a 
larger /3 e ff (Fig. 4a). The heterogeneity in the weighted degrees s in and 
s out is captured by H = cr in a out /(s), where a 2 n and cr 2 ut are the variance 
of the marginal probability density functions P(s in ) and P(s out ) respec- 
tively (Fig. 4b). Finally the symmetry between s m and s out is captured 
by S= ((s in s out ) — (s in )(s out )) /(cr in cr out ), the in-out weighted- degree 
correlation coefficient. This term, — 1 < S < 1 , is positive when nodes 
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Figure 4 | The impact of Ay on resilience. The topological characteristics 
that affect a systems resilience through (3 e s (equation (13)) are: a, the 
network density (s); b, the heterogeneity in degrees or link weights Ti; and 
c, the symmetry 5, capturing the correlations between a nodes in and out 
degrees, d, Phase diagram for mutualistic dynamics in the (s)-TC plane. 

In the resilient phase, the system has a single stable fixed point x H ; in the 
non-resilient phase the undesired x L is also stable. For this dynamics the 
greater /5 e ff is (square size) the deeper the system is in the resilient phase, 
e, The average state of the system (x) versus the average reduction in the 
link weights. For Net 12, the most heterogeneous of the fourteen 
mutualistic networks, we observe an extreme degree of resilience, avoiding 
collapse up to/ w = 97% (blue triangles). A homogeneous network, with the 


same density and TC = 0 loses resilience at/ w = 66% (red circles); hence Ti 
is the source of Netl2’s exceptional resilience. As indicated in d, Netlls 
resilience is rooted in its high density (s). Indeed, we find that for Netl 1, 
both the original (blue triangles) and the homogeneous (red circles) 
networks feature comparable levels of resilience, indicating that Ti has a 
marginal contribution, f, The phase diagram for directed transcription 
regulatory networks is 3D— (s), Ti and S — with the first-order transition 
from a living cell to cell death occurring at ( s ) + STi = 2. The position of 
the S. cerevisiae and the E. coli networks is also shown (triangles), g, For 
both organisms 5< 0, and hence heterogeneity decreases their resilience: 
indeed, the homogeneous networks (red circles) withstand larger 
perturbations than the original networks (blue triangles). 


with large s) n tend also to have a large s° ut . An undirected network 
= Aji is a perfectly symmetric network, as s| n = s° ut , hence S= 1; an 
asymmetric network is one where nodes with a large in-degree tend to 
have a small out-degree, in which case S < 0, tending towards —1 
(Fig. 4c). 

Equation (13) helps us to identify the network characteristics that 
can enhance or weaken a systems resilience. Consider for example the 
resilience of the ecosystem described in equation (5). The structure of 
the resilience function (Fig. 2m) indicates that the greater is (5 e ff, the 
more resilient is the system, enduring larger perturbations before reach- 
ing the bifurcation at /3 c df and risking a transition to the undesired x L . 
Since mutualistic networks are symmetric we have 5 = 1 and 
Ti = cr in cr out /(s) = <t 2 /(s), obtaining /^ e ff = (s) + cr 2 /(s). Hence, equation 
(13) increases both with density (s) and with heterogeneity Ti giving 
rise to a 2D phase space: a resilient phase above the phase boundary 
(s) + Ti = /3 df , and a non-resilient phase below it. In Fig. 4d we show 
all fourteen mutualistic networks and their location on the 
(s)-TC phase diagram, characterizing the source of each networks resil- 
ience. For instance, Netl 1 and Net 12 have comparable /? e ff values and 
hence comparable levels of resilience, indicating that both can with- 
stand comparable levels of perturbation before crossing the bifurcation 
at /3 c e{{ . However, Fig. 4d shows that while the source of Netl Is resil- 
ience is its high density (s), the source of Net 12 s resilience is its high 
heterogeneity Ti. To test this we constructed two homogeneous 
networks with the same densities as Netl 1 and Netl2, but with Ti = 0. 
As expected for Net 12, reducing heterogeneity greatly decreased resil- 
ience (by ~30%), while for Netl 1, whose source of resilience is density 
rather than heterogeneity, eliminating Ti had only a negligible impact 
(Fig. 4d, e). 

In regulatory dynamics the phase diagram has two domains, cor- 
responding to an active phase and cell death. Here resilience increases 


with f3 e ff, as the larger is /3 e ff, the deeper is the system into the active 
state, and farther from the critical transition at /3 d{ (equation (12)). 
Since A y is directed, Ti can both increase or decrease /? e ff, depending 
on the symmetry 5 in equation (13). Consequently, resilience is 
governed by three topological characteristics, where dense, symmet- 
ric and heterogeneous networks are most resilient (large /3 e ff), and 
sparse, antisymmetric and heterogeneous networks are least resilient 
(small /3 e ff). The phase diagram in the ( s)-Ti-S space is shown 
in Fig. 4f, and the transition between states occurs along the 
plane (s) + STi = /3 d{ . For the two regulatory networks we measured 
5= —0.083 (for S. cerevisiae) and S= —0.2464 (for E. coli), both 
negative. Hence here Ti has a negative contribution to resilience and 
a homogeneous network with Ti = 0 would, in fact, be more resilient. 
To test this prediction we compared the resilience function of the 
empirical networks ( E . coli and S. cerevisiae) with that of the equiva- 
lent homogeneous networks in which (s) is preserved and Ti = 0 
(Fig. 4g). Indeed, we find that eliminating heterogeneity increases the 
system’s resilience. 

Complex systems are characterized by an inherently multi- 
dimensional parameter space, giving rise to diverse and unpredictable 
behaviour. Here, by reverting to the natural parameter space (/3- space) 
we exposed the hidden universal patterns of network resilience. The 
origin of this universality is in the separation of the systems dynamics 
and topology. Indeed, in most systems the intrinsic behaviour of the 
components and the nature of the interactions between them are invar- 
iable to perturbations 13 . Perturbations only affect the structure of the 
network, Ay, determining who interacts with whom and how strongly. 
Our formalism reduces Ay into an effective ID system, showing that 
regardless of the specific topology and weights, or the form of pertur- 
bation, the patterns of resilience depend only on the systems intrinsic 
dynamics. The role of the network topology is fully captured by the 
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ID f3 e ff, predicting that density, heterogeneity and symmetry are the 
three key structural factors affecting a systems resilience. They do not 
alter the critical points, but instead push a system far away from these 
critical points, helping the system to sustain large perturbations. This 
separation of structure and dynamics provides us with testable predic- 
tions for the systems response to different perturbations. It also suggests 
potential intervention strategies to avoid the loss of resilience 26-28 , or 
design principles for optimal 29 resilient systems that can successfully 
cope with perturbations 30 . 
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